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ABSTRACT 

We investigate the effect of dust on the observed rotation rate of a stellar bar. The only direct 
way to measure this quantity relies on the Tremaine & Weinberg method which requires that 
the tracer satisfies the continuity equation. Thus it has been applied largely to early-type barred 
galaxies. We show using numerical simulations of barred galaxies that dust attenuation factors 
typically found in these systems change the observed bar pattern speed by 20-40 percent. We 
also address the effect of star formation on the TW method and find that it does not change the 
results significantly. The results presented here suggest that applications of the TW method 
can be extended to include barred galaxies covering the full range of Hubble type. 
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1 INTRODUCTION 

The rate at which a bar rotates, its pattern speed, Q p , is the prin- 
ciple parameter controlling a barred (SB) galaxy's dynamics and 
morphology. Most determinations of this parameter are indirect. 
An often used method is to match hydrodyna mical simulations of 
SB galaxies to observed velocity fields (e.g. IWeiner et al J 1 200 3 ; 
IPerez et all2004l : lRautiainen et all2005t for recent examples). An- 
other commonly used method is to identify morphological features, 
such as rings, with resonance radii to infer the pattern speed. 

The only direct and model-independent technique to measure 
Q p is the Tremaine & Weinberg (1984, hereafter TW) method. 
They show that if a tracer satisfies the continuity equation then 
it is straightforward to derive an expression that relates fl p to 
the luminosity-weighted mean velocity, V, and the luminosity- 
weighted mean position, X, of this tracer. Gas is generally not a 
suitable tracer because it is easily shocked, changes states or is 
converted into stars (e.g. Hernandez et al. 2005; see also Rand & 
Wallin 2004). The assumption underlying the TW method, that the 
observed intensity is proportional to the tracer's density, limits its 
use to relatively dust and gas-free systems, i.e. early-type galax- 
ies. The model-dependent techniques, on the other hand, rely on 
the presence of gas and are therefore restricted to late-type barred 
galaxies. There is therefore no overlapping range in Hubble type 
where both methods have been applied and compared. 

The number of successful applications o f the TW 
method was initially rath er limited: N GC 936 dKentl 1 19871 : 



iMerrifield & Kuiikerl 1 19951) NG C 4596 dGerssen etalTll999.. 
NGC 1023 dDebattista et alj 20021) . More recently, larger samples 
have been obtained: Aguerri et alj d2002l) present a study of five 



SB0 galaxies and Gersse net alj d2003h study four early-type 
barred galaxies. The first successful application of the method 
using integral fie ld (Fabry-Perot) absorption-l ine spectroscopy 
was presented by iDebattista & Williams! d2004l) for NGC 7079, 
while the first direct detection of kinematic decoupling in a 
double-barred g alaxy, NGC 295 0, using the TW method was 
presented by ICorsini et alj (2003]). A surprising result emerging 
from the se studies is that bars r otate as fast as they physically can 
(see also lRautiainen et alj2005h . 

Studies of dust extinction in disc galaxies have a long 
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and checkered history (e.g. iHolwerda et al]|2005L and references 
therein). In this paper we seek to quantify the effect of dust absorp- 
tion on the TW method. To that end we will u se th e A^-body simu- 
lations of barred galaxies by Debattiste] d2003l) and lDebattista et all 
( 2006) and implement dust extinction by adjusting the weights of 
individual particles. Our motivation for this study is to explore how 
important dust is for TW measurements of both early and espe- 
cially late-type galaxies. In the latter, star formation can be a poten- 
tial problem as it violates the basic assumption of the TW method. 
We therefore also explore the ability to measure fi p in a hydrody- 
namical simulation of a barred galaxy that includes star formation 
dDebattista et alj2006l) 

Barred galaxies often display prominent dust lanes along the 
leading edges of the bar. Numerical simulations associate the dust 
lanes with shocks in the gas streaming along the length of the bar 
( Atha nassoulal 1 1 99*21) . Observations of velocity jumps across dust 
lanes, as well as the detection of dust lanes in radio continuum, 
confirm this interpretation (e.g. lMundell & Shondll999 ). With the 
introduction of efficient spectrometers operating in the IR the possi- 
bility now exists to carry out TW measurements in late-type galax- 
ies, eliminating an important observational bias. We therefore test 
whether such measurements would be compromised by dust extinc- 
tion and star formation. 
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Figure 1. left: panel schematically illustrates the straight dust lane geometry (cf Figure 1 in lAthanassouk]jl992l) ') used in our numerical experiments. The bar 
is represented by the dashed ellipse in this cartoon, right: This cartoon shows that TW measurements of the bar at a positive PA (solid lines) are not equivalent 
to the same observations at a negative PA (dashed lines). This is because the long-slits (hatched areas) intersect the dust lanes, represented by the straight line 
segments, at different locations depending on the sign of the bar position angle, which leads to a change in sign of the error in Q. 



Our results show that dust lanes on the leading edges of a bar 
change the TW derived value of fl p depending on the geometry 
of the observations. This can lead to measured pattern speeds that 
differ by 40 percent or less from their true value for realistic dust 
lane models. As this change can be rather small and because the 
extinction can be observationally estimated, it suggests that the TW 
method can be applied to the full range of Hubble type. 

The paper is organized as follows. We summarize the TW 
method and describe the models to which we apply it in Sections 2 
and 3 respectively. In Section 4 we describe our models for the dust 
lane and the diffuse dust disc. The model predictions are presented 
in Section 5 and in Section 6 we compare the models to data and 
discuss the implications of our modeling. 



2 TW METHOD 

For any tracer population that satisfies the condition of continuity, 
it can be shown that 

V = XQ p smi. (1) 

Here, X = J h(Y)XE dXdY and V = / h(Y)V los T, dXdY are 
luminosity-weighted positions and velocities respectively. The sur- 
face brightness of the tracer population is E, and V\ os is the line-of- 
sight velocity minus the systemic velocity. The galaxy's inclination 
is given by i and coordinates along its major and minor (disc) axes 
by (X, Y). The weighting function, h(Y), can be chosen arbitrarily 
(TW). For several slit observations (i.e. h(Y) oc 8(Y — Y a set )), 
plotting V versus X produces a straight line with slope Q p sin i. 
In our experiments, we will refer to the slopes obtained from fit- 
ting such straight lines as Q sin i, in order to distinguish from the 
intrinsic pattern speed f2 p obtained from the time evolution. 

We apply the TW method to the TV-body models of barred 
galaxies described in Section|3]as done by Debattista (2003, here- 
after D03). In our measurements, except where noted, we set the 
inclination i = 45° and the projected angle between the bar major 
axis and disc major axis, PAb ar = 45°, implying an intrinsic angle 
of 54.7°, which is a favorable orientation for the TW method. 



As in D03 we use 11 slits each with a width of 0.1 (ini- 
tial) disk scale-length s (the units of the simulation are described in 
Debatti sta et al . 2006) to derive Q. Tests with slits half and twice as 
wide return nearly identical results. Applications of the TW method 
to real data generally use between three and five long-slit spectra. 
(Integral field observations could soon change this.) A typical ob- 
servational setup places one slit on the major axis and two slits on 
opposite sides either near the bar edges, or at half the bar radius. 
Tests with such configurations in the simulations did not yield val- 
ues of Q that were significantly different from the results derived 
with all 1 1 slits. The strongest constraints on the slope come from 
slits located about half a bar radius from the galaxy center. The 
leverage on the slope of the fit derived from spectra beyond this 
point is decreasing with radius, as is their signal-to-noise ratio (cf. 
D03). Progressively removing pairs of the outermost slits from the 
set of 11 and repeating the measurements has little effect on Q, 
although it increases the scatter. 



3 BARRED GALAXY MODELS 

In their study of the effec t dust has on the observed kinematics 
of discs, iBaes et a 1. (2003) used analytical distribution functions 
for the intrinsic stellar kinematics. However, such methods are ill- 
suited for modeling barred galaxies. In this paper we explore the 
effect of dust on barred systems using TV-body and TV-body+SPH 
models. 

3.1 TV-body models 

Our four TV-body si mulations TW1-TW4 were evolved on a polar 
cylindrical grid code dSellwood & Valluril 1997b . The radial spacing 
of grid cells in this code increases logarithmically from the cen- 
ter and reaches to well past the edge of the model discs. We used 
Fourier terms up to m — 8 in the potential, which was softened 
with the standard Plummer kernel. Time integration was performed 
with a leapfrog integrator using a fixed time-step. 

These models use rigid halos; this allowed us to sample the 
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disc (and bulge, if present) component at high mass resolution. 
A full descript i on of these models can be found in D03 and 
iDebattista et al.l d2006t) . Briefly, the rigid halos were represented 
either by a spherical logarithmic potential or by a Hernquist halo. 
In all models, disc kinematics were set up using the epicyclic ap- 
proximation to give a constant Toomre Q. We assume that all the 
late-type models are bulgeless and only include a bulge component 
(using the method of Prendergast & Tomer, 1970) in the early-type 
model presented in D03. We use units in which the initial radial 
scale length of the disc particle distribution is 1; in models TW1- 
TW4 the vertical scale height is then 0.1 in these units. Except 
where noted, we report results using these natural units. The result- 
ing axisymmetric systems all formed rapidly rotating bars. 

The model TW5 is an unpublished simulation with an expo- 
nential disk + dominant live dark matter halo evolved with PKD- 
GRAV. Because the halo was live and dominant in this simulation 
the bar ends well inside its co-rotation radius, i.e. it is a slow bar 
(Debattista & Sellwood 1998, 2000). This also causes the bar to be 
larger and stronger than any of the others. 



3.2 7V-body+SPH 

In late type galaxies the rate of star formation (hereafter SF) can 
be substantial. Thus stars violate the required condition of (source- 
free) continuity. To quantify the influence of SF we also include an 
7V-body+SPH model of a barred galaxy in ou r analysis. This model, 
described in detail in lDebattista et alj (2006), included a live halo 
and was evolved us ing a tree cod e (GASOLINE). The star formation 
recipe is the one o f lKatzl dl992l) . In our model the SF is most pro- 
nounced in the central regions of the system and in the bar itself, 
consistent with numerous observational results, (e.g. Sheth et al. 
2002). 

We measured Q, p directly from the time-evolution; since the 
dark matter halos in these models are rigid, f2 p remained constant 
after an early phase of evolution. Table [TJlists the properties of the 
bars, including the pattern speeds fi p , of all the models. 



3.3 Individual models 

The models were selected to have a variety of morphologies and 
pattern speeds. We generally selected late times in the simulations 
to ensure that the models are well settled. Here we use model units. 

TW1. This model of an early- type galaxy was described in 
D03, who used it to study the effect of position angle errors on 
TW measurements. We use the model at t = 2.48 Gyr, by which 
time a strong and stable bar had formed. Only weak spirals are 
present in this model, resulting in a minimal interference with the 
determination of Q p . This model has 4 x 10 6 equal-mass particles 
in bulge+disc. 

TW2. The evo lution of this model is discussed in 
IDebattista et al.l ([2006), where it is referred to as model L5. We use 
model TW2 at t — 7.44 Gyr at which time its spiral is relatively 
weak. It is a late-type model with 7.5 x 1 6 particles. 



TW3. Model TW3 is described in Debatti sta et alj i 2006). 
where it is referred to as model H2, and in lDebattista et alj (2005), 
where it is referred to as R7. We use t = 2.48 Gyr at which time 
spirals are rather weak. This late-type model has 4 x 10 6 particles. 

T W4. This late-type model is described in Debattista et al. 
d2006h where it is referred to as model S 1 . We use model TW5 at 
an early time, t = 1.24 Gyr, when the spiral arms were prominent. 
This model contains 7.5 x 10 6 particles. 



TW5. Slow bar model, with 0.75 x 10 particles in the disk 
and 2 x 10 6 in the halo. We use this model at t ~ 10.8 Gyr. Unlike 
the other models, TW5 is not yet published. 

TW6. This mode l includes star formation. It is described in 
IDebattista etafl (2006) where it is referred to as NG3. We use this 



model at t 
particles. 



3.48 Gyr. At this time the model contains 0.7 x 10 



4 DUST MODELS 

Observations of dust in disc galaxies show that the dust dis- 
tribution in these sys tems follows a double exponential profile 
(Wain scoat et alj 1989 ) and it is conventional to m odel them as such 
(e.g. lBaes et alJl2003l : IValotto & Giovanellill 20041) . Thus we model 
diffuse dust discs as: 



D disc (R, z) = D ,disc e 



-R/h r 



MA*,, 



(2) 



where /iodise and /i z ,disc are the radial scale-length and the vertical 
scale-height of the dust disc respectively. The central dust density 
is given by Do.disc 



4.1 Dust Lanes 

Late-type barred galaxies often exhibit prominent dust lanes run- 
ning from th e end of the bar toward the center of the galaxy. 
Athanassoufa] j 1992b notes that these dust lanes generally fall in 
two categories: (i) straight dust lanes which are parallel to each 
other but make an angle with the bar major axis, and (ii) curved 
dust lanes which are parallel to the bar major axis but offset toward 
the leading edges and with their concave side facing the bar major 
axis. We implemented both morphologies in our experiments. Since 
the results obtained with the curved dust lane morphology and the 
straight dust lane morphology are qualitatively and quantitatively 
similar we concentrate only on the latter in this paper. 

We implement dust lanes by defining rectangular areas on the 
bar's leading edges that mimic Athanassoula's straight dust lanes. 
The dust lane geometry is fully specified by four parameters: the 
width of the dust lane (tOdi), the angle the dust lane makes with the 
bar major axis (sai), and the minimum and maximum extent of the 
dust lane (x m i n and x max respectively) along the x (bar) axis. In 
practice we fixed Kmax to the radius of the bar and the dust lanes 
start on the bar's major axis. A cartoon representation of these var- 
ious parameters is shown in the left panel of Figure [T] Note that 
x m in can be negative, i.e. the dust lanes can extend past the center 
of the bar. In this case, the inner parts of th e dust lanes are some - 
times observed to wind around the bulge, see lAthanassoulal dl992l) . 
The sign of PAbar is important, as illustrated in the right panel of 
Figure Q] because the slits used to measure the TW integrals inter- 
sect the dust lanes at different locations for positive and negative 

PAbar • 

We assume that the dust distribution within the dust lanes can 
also be described by a double exponential model: 



Aane(#, z) 



£>0,lane e 







if in dust lane, 
otherwise. 

(3) 

Although the dust lane model has a similar double exponential pro- 
file as the dust disc model, the scale parameters of the two compo- 
nents will in general be different. 
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Figure 2. Examples of the projected distribution of particle weights in our numerical models shown on a logarithmic stretch. The models are shown here 
face-on to highlight their morphological features. Dust lanes have been included to illustrate how dust attentuation is implemented in our simulations. For 
clarity the dust lanes have an exaggerated extinction (A v = 15) and width. 
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Table 1. Our numerical models have bar semi-major axis ag The values of ag were estimated by eye. We measured the pattern speed, f2 p , from the time 
evolution. To obtain the dust-free TW method values, Slo, we orientated the models to i = 45° and projected PAb al . = 45°. The difference between Q p and 
Qq is listed under column Error. 



Run 


as 




n 


Error 


Reference 




kpc 


km s~ 1 kpc -1 


km s kpc -1 






TW1 


4.5 


23.3 ±0.8 


26.5 ± 0.8 


±14% 


Debattista (20031 


TW2 


7.8 


15.0 ±0.2 


15.9 ±0.1 


±6% 


Debattista et al. (2006) 


TW3 


3.8 


22.8 ±0.2 


24.2 ±0.2 


±6% 


Debattista et al. (2006) 


TW4 


3.3 


34.6 ± 1.1 


38.0 ±4.7 


±10% 


Debattista et al. (2006) 


TW5 


8.0 


8.8 ±0.4 


9.5 ±0.5 


±8% 


unpublished 


TW6 


5.9 


26.4 ±0.2 


25.3 ± 1.1 


-4% 


Debattista et al. (2006) 



4.2 Particle weights 

In a particle implementation of the TW method, the integrals in 
Eqn. Q] are replaced by sums over particles. Each particle has its 
own weight, w. In the collisionless iV-body models (TW1 to TW5) 
all particles are coeval and have equal mass. We therefore assign 
the same intrinsic (i.e. before applying dust attenuation) weights to 
all particles in these models 

In model TW6, however, the particle luminosities, and hence 
their, weights will depend on their age. To account for this we scale 
the intrinsic weights of all particles by their M/L ratio in this model. 
We do not a priori know the luminosity of a particle. We only know 
its mass M p and time of creation t p . To calculate the luminosity of 
a particle we follow the procedure first described by Tinsley (1973, 
see also van der Kruit 1989) and summarized in the appendix. 

When dust is present, particles obscured by dust contribute 
to the integrals with lower weights. Unlike a foreground screen 
model, the dust extinction in our models varies with position within 
the disc. For each particle we calculate the amount of intervening 
dust by integrating the projected dust distribution along the line- 
of-sight. The weight for each particle is then given by Wi — e~ Ti , 
where the optical depth n is obtained from 




assuming a unit mass absorption coefficient. This integral is evalu- 
ated using a Gauss-Legendre 64-point quadrature formula from the 
FORTRAN NAG libraries. To guard against numerically unstable 
solutions we subdivided the integration interval into ten subinter- 
vals to improve the accuracy. Examples of the weight distribution 
in each of our six models are shown in Fig. [2] They are shown here 
face-on to better illustrate their morphological features. 

Following convention, we specify the amount of extinction in 
our models as the face-on optical depth for a particle infinitely far 
behind the disc, to. For the dust disc models, we place this particle 
behind the center of the disc, while for the dust lanes we place it 
at a point halfway along the dust lane. Equivalently, the extinction 
can be expressed in Ay magnitudes 

Ay = 2.51og 10 er = 1.086r . (5) 
4.3 Code test: A dust disc 

IValotto & GiovaneliU |2004, hereafter VG04) presented an analyti- 
cal model of a diffuse dust disc and used this model to demonstrate 
that the presence of dust lowers the slope of the observed inner 
rotation curve at large inclinations. As a test of our dust extinc- 
tion code we reproduce their result using a simple axisymmetric 
double-exponential model of a cold rotating disc with Rd = 1.0 
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Figure 3. To test our dust code we reproduce the analytical results of VG04 
and derive the major axis rotation curves in an axisymmetric particle data 
set. We added a diffuse dust disk with a face-on central optical depth of 
ro = 3.0 and derive the rotation curve at various inclinations. The projected 
amount of dust extinction varies with inclination leading to more slowly 
rising rotation curves. From top to bottom the inclination of the rotation 
curves varies between 70° and 86°. At the largest inclination, the rotation 
curve does not reach its asymptotic values before the edge of the A^body 
data (cf Figure la of VG04). 

and Zd = 0.1 respectively. The circular velocities are implemented 
using equation 2 of VG04 and their 'luminosity class 1' parame- 
ters. We sample this model randomly by 10 6 particles. Dust in this 
model is described by a diffuse disc with geometric parameters that 
are identical to those for the particle distribution: fo^disc ~ Rd 
and /i 2 ,disc ~ Zd- In order to facilitate a direct comparison with 
the results of VG04 we use the same dust extinction in the disc 
(ro = 3.0) as they do in their Figure la. 

We measured the rotation curves in our model by radially bin- 
ning the particle velocities in a slit along the major axis. In each ra- 
dial bin the line-of-sight velocity distribution (LOSVD) of the par- 
ticles has a sharp cut-off at high velocities and a tail toward lower 
velocities. We take the cut-off velocity as the measure of the cir- 
cular velocity in each bin. The resulting rotation curves are shown 
in Fig. [3] As these curves are derived from a finite number of par- 
ticles they are not as smooth as the analytical results of VG04 but 
their qualitative and quantitative behaviour is quite similar: inter- 
nal extinction decreases the slope of the inner rotation curves. This 
agreement with the results of VG04 strengthens our confidence in 
our extinction code. 



5 RESULTS 

The TW-measured value of the bar pattern speed, Q, can be af- 
fected not just by dust but also by spiral structure. We quantify the 
effect of dust using the ratio f2/f2o, where Qo is the TW measured 
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Figure 4. Illustrations of the effect of varying the dust lane model on Q. The contours show the ratio Q/Qq. In each panel we evaluate this ratio at a fixed value 
of the extinction (Ay ~ 8) and with a dust lane slope of s c \i = —0.3. Unlike the axisymmetric dust disc model, the sign of PAb ar matters (see Figure[TJ. The 
plotted contours are for PAb ar = +45° (solid) and PAb ar = —45° (dashed). left: The ratio Q/CIq for varying dust lane width (uidl) an d minimum dust lane 
radius (rr m j n ). centre: The behaviour of the ratio Q/Qo as a function of viewing parameters. In this panel the dust lane model parameters are fixed at x m ; n = 
—0.1 and M>dl = 0-1. right: The dependence of the ratio on the scale length parameters of the dust lane model. For reference, the initial scale lengths of the 
disc particle distribution are R c \ = 1.0 and z c \ = 0.1. The cross in the left and middle panels marks the approximate location of NGC 4123 in these panel 
(see Section|6]but note that the extinction in NGC 4123 is only Ay ~ 3). 



value in the absence of dust. To quantify the effect of spirals, we 
compare Qq with S7 P , the pattern speed inferred from the time evo- 
lution of the models, in the absence of dust. The resulting values of 
S7 for each TV-body model are listed in Table[TJand compared with 
S7 P . The TW measurements generally differ from the time evolution 
values by 15% (as already noted by D03). Model TW5 included 
strong spirals in the disc; despite this, S7 differs from S7 P by less 
than 10 percent. 

We conclude that in the presence of a strong bar, spirals only 
perturb TW measurements rather weakly. This is presumably be- 
cause spirals are at large radius and therefore in lower density parts 
of the disc. Moreover, because spirals wind around the galaxy cen- 
tre, their contribution to any slit partly cancels out. Therefore, spi- 
rals are not expected to present a serious problem for TW measure- 
ments in late-type galaxies with strong bars. (But note that, other 
than obscuration, we give all particles equal weight; in a real spi- 
ral star formation would increase its relative contribution in blue 
bands.) In the remainder of this paper, we will compare the mea- 
surements with and without dust directly, rather than with the time 
evolution value. 



5.1 Diffuse dust disc 

We first briefly consider the effect on SI of a diffuse dust disc. 
We examined the ratio S!/S7o resulting from varying the parame- 
ter Do, disc, and therefore Ay. Typically observe d values for Ay 
range from 0.5 to 4.0 mag (Holwerda etai]|2005h . The change in 
S7 is only five percent for Ay ~ 3. Even at Ay = 8, the change 
is less than 15 percent. We also examined the dependence on the 
assumed /i_R.disc and ft z ,disc of the dust disc at a constant to = 5. 
Unsurprisingly, the value of S7/ flo depends only weakly on either 
scale parameter at fixed to once they become of order of, or larger 
than, the scale of the stellar system itself. Even for to = 5, the error 
in Q does not exceed 10 percent. 

In real galaxies, dust within the bar radius is likely to be swept 
into a dust lane. To simulate this we tested what the effect of a dif- 
fuse d ust disc with a hole within the bar radius is (cf Marshall et al. 
2006). We set h R , d i S c = 10, /i z ,disc = 0.1, and r = 10. As in- 



tuitively obvious, such a dust disc has a negligible (< 2 percent) 
effect on SX 

Diffuse dust discs lead to a small errors in Q, of order 10 per- 
cent. Existing TW measurement of early-type barred galaxies that 
lack dust lanes but possibly have diffuse dust discs are therefore un- 
likely to have been affected substantially by dust. All bars with TW 
measurements so far are consistent with being fast, in the sense that 
they end near their corotation radius. This fi nding supports low dark 
matter densities at the centers of g alaxies ijDebattista & Selfwoocj 
ll998l ; ISellwood & Debattistdl2006T) . 



5.2 Dust lanes 

We therefore ignore the dust disc contribution in our dust lane mod- 
els and start by exploring the dust lane geometry and the dust pa- 
rameters using model TW1. This is a strongly barred model. As 
strongly barred galaxies are observ ed to have flat surface bright- 
ness profiles along their major axis dElmegreen et al J 19961) . we set 
/iji.iane = 10 to obtain a nearly constant radial dust lane distribu- 
tion. For the vertical scale height of the dust lane we set a value of 
foz.iano = 0.1, which is comparable to the initial stellar distribution. 

Using this distribution for the dust lane and Do.ianc = 100 
(corresponding to an extinction of Ay ~ 8) we evaluate the ratio 
S7/S7o for varying s<u, Wdi and a; m i n - The solid lines in the left 
panel of Figure [4] illustrate the variation in this ratio as a function 
of Wdi and x m i n for constant Sdi = —0.3. For all combinations of 
parameters the derived ratio is larger than one. Dust lanes therefore 
lead the TW method to overestimate the pattern speed by 10 to 40 
percent. Variations in the width have the largest impact on Q, as they 
affect the geometric area of a dust lane the most. 

As noted in Section |4~T1 the (observationally known) sign of 
PAbar is important for the dust lane model. We illustrate this in the 
same panel using dashed contours for PAb ar = —45°. The effect 
of a negative PA is that S7 becomes smaller than Qq rather than 
larger as it does for positive PA. Dust lanes can thus also lead the 
TW method to underestimate the pattern speed. For realistic values 
of the model parameters (e.g. the NGC 4123 values marked by the 
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cross, see also section^ the percentage change in Q/ f2n is roughly 
independent of the sign of PAbar. 

We repeated the exercise for a range of dust lane slopes brack- 
eted by —0.4 ^ sai ^ —0.1. The behavior of fi/fi u is qualitatively 
and quantitatively similar for all but the most shallow slope. When 
s<H = —0.1, the errors in Q, are rather small. This happens because 
the dust lane is no longer near the leading edge of the bar but closer 
to the bar major axis. However s<ji= —0.1 bears little resemblance 
to realistic dust lanes. For most combinations of the dust lane ge- 
ometry parameters the observed change is Ss 25 percent. 

5.2.1 Viewing parameters 

Inclinations of ~ 45° and projected angular separations between 
the bar major axis and the disc major axis of ~ 45°, such as we 
used thus far, are favorable orientations for TW measurements. In 
the central panel of Figure [4] we examine the variation in n/fio as 
a function of inclination and bar position angle. (The dust lane pa- 
rameters used in this calculation are s m i n = —0.1, Wdi — 0.1 rela- 
tive to the bar semi major axis and sai = —0.3, which correspond 
to the approximate values observed in the SBc galaxy NGC 4123 
scaled to our model.) The error in fi reaches a maximum when 
the bar is close to the minor axis at extremely low or high incli- 
nations. Galaxies with such orientations are not usually selected 
for TW measurements: at large inclinations the projected disc on 
the sky is too narrow to employ the multiple parallel long-slits re- 
quired by the TW method while at low inclinations the observed 
velocities, Vi nt sin i, will generally be too small to measure reli- 
ably. Around i = 45° and PAbar = 45°, the results are indepen- 
dent of bar position angle and depend only weakly on inclination. 
For 30° < i < 60° and 20° < PA bar < 80°, these errors are less 
than 25%. 

5.2.2 Scale length parameters 

The observed pattern speed depends only weakly on variations in 
ftfijanc- As our initial /i_R,iane is ten times that of the particle distri- 
bution, even large variations around this initial value do not lead to 
significant changes in Q. The dependence on /i 2! i ane is somewhat 
stronger as illustrated in the right panel of figure [4] Once /iK,i an e 
or /i z ,i an e become larger than the corresponding quantities for the 
stellar distribution, fi/flo becomes independent of these parame- 
ters. 

5.3 Different TV-body models 

In the previous sections we explored the behavior of Q, on the dust 
model parameters with only model TW1, largely to facilitate com- 
parisons with D03, who used the same model. We carried out a 
similar parameter space study for all the other TV-body models pre- 
sented in Section [3] The same dust lane geometry (with Xmin and 
Wdi rescaled to ob) and dust distribution used in Section [5~2l was 
used on each TV-body model. Only the maximum radii of the dust 
lanes (x max ) differ to account for the different bar lengths in each 
model (see Table [TJ. The dust lane widths used here are fairly nar- 
row, about five percent of the bar semi-major axis length, to re- 
semble the shapes of observed dust lanes. These yield qualitatively 
similar results. A quantitative comparison of the different TV-body 
models is shown in Figure|5]where the ratio fl/fio as a function of 
extinction is plotted for each model. The main difference is in the 
errors in model TW5 at positive P Ab ar '■ the error for this slow bar 
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Figure 5. The ratio Q/Qo of the bar pattern speed observed with the TW 
method to the intrinsic bar pattern speed as a function of dust lane extinction 
Av (defined face-on as described in i|4.2t . The different curves show the 
behaviour of this ratio for the different TV-body models. Two curves are 
shown for each model. Curves with a ratio > 1 are at PAb ar = +45° , while 
curves with a ratio < 1 are at PAb ar = —45°. The dust lane geometry and 
dust distribution are the same for every model: w^i/ag = 0.1, Sdl = —0.3, 
x m i n /aB = —0.1, hji i anc = 10 and /i 2 j ano = 0.1. For clarity, the lc 
eiTors are shown for model TW1 only. 

is a factor two larger than in the other models. However at negative 
PAb ar the error in this model is not much different than the other 
models. 



5.4 Star formation 

The models that we explored so far do not take star formation into 
account. However, in late type systems the SF rate can be large 
enough that departures from continuity become too large. Bright 
young stars will contribute disproportionately to the TW integrals 
X and V (see section [2]l and hence may significantly change the 
measured pattern speed. To quantify the effect of star formation on 
the TW method we use model TW6. The result is shown in Fig. [6] 
and is qualitatively similar to the results derived from model TW1 
to TW5. Quantitatively, the impact of dust in TW6 lies between 
the values in models TW1-TW4 and TW5, with an value that 
differs from its true values (at A v — 3) by 15 percent. Scaling the 
intrinsic weights of the particles by their M/L ratio (section [3~2t 
is a crucial part of model TW6. As a test we therefore derive the 
pattern speeds in this model without any such scaling (dotted line 
in Fig.(6]l. Without this scaling the effects of dust are comparable to 
the slow bar model TW5 with f2 deviating from the intrinsic value 
by 21 percent. 

To examine this more closely we divided the particles in two 
groups, those particles that were present in the simulation from 
t = ('old') and the ones that were created subsequently by star 
formation ('new'). The pattern speed derived from younger parti- 
cles appears to be less affected by dust. In fact, subdivisions into 
younger ages show even shallower curves. Interestingly, the dust 
free values of fio derived for the different age bins vary by less 
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Figure 6. The ratio f2/f2o as a function of dust extinction for the model 
that includes star formation (TW6). The result is not very different from 
the models shown in Fig [5] For comparison model TW1 is overplotted as 
the solid line with error bars. Without M/L scaling (see text) dust affects 
model TW6 significantly more than our previous models. The other curves 
show the behaviour of Q in TW6 for particle subsets in age. The younger 
the particles the less the effect of dust appears to be. 

than eight percent from the full model. This may be somewhat for- 
tuitous as there are fewer particles at younger ages to derive the 
TW value. Their spatial extent is also smaller as these particles are 
mostly confined to the bar region and the nucleus. 

This behaviour is probably due to the similar location of dust 
lanes and SF regions. If star formation occurs or is enh anced on the 
leadin g edges of a bar, as is frequently observed, e.g. She th et all 
j2002l) , then absorption by dust may be partially mitigated by the 
enhanced luminosity of the young particles. SF is indeed observed 
along the bar in this model. 



6 DISCUSSION 

The experiments presented in this paper show that dust lanes on 
the leading edges of a bar tend to increase the TW derived value 
of the bar pattern speed, Q, when PAbar > and decrease it when 
PAbar < 0. The changes in fl are relatively modest for realistic 
dust lanes (i.e. Av 55 3), between 8-25 percent (cf. Figs [5] and [6) 
and increasing to 20-40 percent at unrealistically large Av ~ 8 
or if the bar is slow. Somewhat surprisingly, the simulation that 
includes star formation (TW6) does not change this conclusion in a 
significant way. Dust lanes alone therefore pose no serious problem 
when trying to extend the TW method to late-type barred galaxies. 

Figure [7}i illustrates why the measured Q changes when dust 
lanes are present. Shown by solid lines in this figure are the spatial 
and velocity profiles along three slits in model TW1 in the absence 
of dust. The dashed lines show the same profiles but for a model 
that includes leading edge dust lanes when PAb ar = +45°, while 
the dotted lines are for PAbar = —45°. The spatial profiles clearly 
show the asymmetric behaviour expected for this configuration (cf 
Figure[Tp). The measured pattern speed Q, follows from the slope of 
the linear fit to the means of the velocity and spatial profiles. The 



dust distribution affects both the spatial and the velocity profiles. 
If dust would alter only the spatial profiles but leave the velocity 
profiles unchanged then the slope of the linear fit would flatten, 
leading to slower pattern speed when PAbar = +45°. However, 
the associated change in V due to the dust more than offsets this 
effect and the bar pattern speed therefore appears larger with dust 
than without. 

The bars in models TW1 to TW4 are classified as 'fast'. That 
is to say that the ratio of the co-rotation radius to the bar ra- 
dius in these models is close to but greater than 1.0. This ratio 
cannot be smaller than 1.0 as that woul d violate self consistency 
(Contopoulos & Papayannopoulos 1980). Applications of the TW 
method to real data for early-type galaxies have not revealed slow 
(i.e. with a ratio > 2) bars. We nevertheless included a slow bar, 
model TW5, to examine the effect of dust in such a case. In this 
model the impact of dust on fl is larger than in the fast bars (Fig.[5} 
This is not surprising as the absolute change in Q between succes- 
sive steps in A v is, to within a factor two, independent of the details 
of a particular model. In other words, the factor by which dust acts 
to change Q is similar for all models and hence makes more of an 
impact on a slow bar. A similar behaviour is observed in model 
TW6 (Fig. [6j. The ratio of bar radius to corotation radius is about 
1.5 indicating that this bar is slow-ish. The observed pattern speed 
is indeed more affected by dust than in models TW1-TW4. 

An important concern with TW measurements i n late-type 
system s is the requirement of a well constrained PA. iDebattistal 
d2003l) demonstrates that a PA error o f 5° can result in pattern spee d 
errors by up to 100 percent (see also lDebattista & Williamsll2004l) . 
In real galaxies, non-axisymmetric features can lead to uncertain- 
ties in th e derived position angles (a nd inclinations) that are on av- 
erage 5° ( iBarnes & Sellwoodll2003l) . It is therefore crucial to select 
systems with well constrained PAs when applying the TW method 
to real late-type galaxies. 

IWeiner et al.1 ( 1200 II) presented a B — I color map of the SBc 
galaxy NGC 4123 that clearly delineates the dust lanes on the lead- 
ing edges of its bar. They found t hat (B — I) = 2.3 in th e dust 
lanes. With the extinction law of iRieke & Lebo f sky (198^) this 
colour implies an average extinction in the dust lanes of Av = 2.7 
magnitude. We checked this value using an HST archival WFPC2 
F450W band image of NGC 4 1 23 (SNAP-9042 PI: Smartt) to com- 
pare a surface brightness profile across the dust lane (i.e. parallel to 
the bar minor axis) in the HST image to a similarly placed slit in 
model TW1. We us e Av = 3 to match the value implied by the 
IWeineretalJfeOOll) colour map. As the exposure time of the HST 
image is only 160s we increase the S/N of the surface brightness 
cut by summing over 100 pixels along the dust lane. This corre- 
sponds to a length of about 10 arcsec over the dust lane or a fifth of 
the bar radius (Weiner et al. 2001). The ragged line in Figure [7}) 
shows the resulting surface brightness profile. The dust lane itself 
is clearly noticeable as the depression near y ~ — 5 arcsec. The 
smooth bold line represents a similar cut along the model image 
averaged over a similar fraction of the dust lane. The model profile 
is scaled ad hoc to match the peak in the NGC 4123 data: the re- 
sulting profile fits the observed one remarkably well. Using this ex- 
tinction, the models shown in Figures[5]& 6 predict a ~ 15 percent 
change in the observed bar pattern speed. If, on the other hand the 
observations were carried out in /("-band using the CO bandhead, 
the extinction of Ay = 2.7 magnitudes translates to Ak = 0.112 
(e.g. lRieke&Lebofekvlll985l . their table 3). This extinction would 
lead to insignificant errors in f2 from obscuration. 

Our results suggest that the application of the TW method can 
be extended to later-type barred galaxies. We therefore propose that 
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at this time a test case of a TW measurement in a late-type galaxy is 
warranted. A useful approach would be to try a TW measurement 
on a bar with pattern speed measured by hydrodynamical modeling 
dWeiner et aljhoOlHPerez et alj|2004) . 
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APPENDIX 

We know the mass M p and time of formation t p of all particles 
in model TW6. To calculate the luminosity of a particle we follow 
Tinsely (1973). We begin by assuming that the stars that make up 
a particle (typical mass ~ 1O 5 M ) are distributed according to a 
Salpeter IMF: f(M) = CM~ 2 35 . The number of stars in a par- 
ticle, C, can be found by setting Af to txC= M p . The total mass 
is 

Aftot = / f(M)dM 

J Mia 

where A/ !o = 0.1 M and M up = 100 Af . We can now calculate 
the total luminosity of this particle by assuming that the luminosity 
along the main sequence scales as M a where a — —4.1. Then, 

rM up 

Ltot =C M a M~ 2 ' 35 dM 

at the moment the particle is created. 

At later times the more massive stars will have turned off the 
main sequence. The upper limit for a star on the main sequence is 

t 1/7 

given by M up = 10 Q yr where 7 = —2.5. By using this value 
for M up in the L tot equation the luminosity of a particle at any time 
t can be derived. Stars with masses larger than M up go through a 
luminous but short lived giant phase. We in clude the contributio n 
of giants in our calculation of Ltot following lvan der Kruitl jl989h . 
But only in the older particles do they make a contribution compa- 
rable to the main sequence stars. 
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NGC 4123 HST - F450W 




Figure 7. left: Cartoon illustrating the effect of leading edge dust lanes on the derived pattern speed. The solid lines show the results derived from the TV-body 
model using three slits in the absence of dust. The dashed lines show the results obtained when leading edge dust lanes are included in a model with PAt, ar = 
+45°. The dotted lines show results derived using the same model but with PAt, ar = —45°. The central panel plots the mean line of sight velocities, (V) or V, 
versus the mean position, (X) or X, for each slit. The slope of the linear fit gives Q sin i. The dust lane modifications to the velocity and luminosity profiles 
(shown in the sub panels) lead to either a steeper or a shallower slope of the fit depending on PAt, ar (see Figure[TJ. The velocity profiles from top to bottom 
and the luminosity profiles from left to right are derived from slits located at Y = —0.4, Y = 0.0, and Y = 0.4, where Y is defined as in Figure[T] right: 
The ragged line shows a cut across the dust lane in an HST-WFPC2 image of NGC 4123 in the F450W filter. To improve the statistics the profile is summed 
over 100 pixels (~ 10 arcsec or about a fifth of the bar radius) along the dust lane. The dust lane itself is clearly noticeable as the depression at y as —5 arcsec. 
The smooth bold line shows a similar cut in model TW1 at a dust lane extinction of Ay = 3.0 and summed over approximately the same fraction of the dust 
lane. The model profile is scaled in an ad hoc fashion to match the peak of the observed profile. The agreement with the observed profile is remarkably good. 
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